Unevenness of Loop Location in Complex Networks 
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The loop structure plays an important role in many aspects of complex networks and attracts 
much attention. Among the previous works, Bianconi et al find that real networks often have fewer 
short loops as compared to random models. In this paper, we focus on the uneven location of loops 
which makes some parts of the network rich while some other parts sparse in loops. We propose a 
node removing process to analyze the unevenness and find rich loop cores can exist in many real 
networks such as neural networks and food web networks. Finally, an index is presented to quantify 
the unevenness of loop location in complex networks. 
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Introduction. In the last decade, the researches on 
complex networks have rapidly developed. At the same 
time, the loop structure has attracted much attention. 
Loops are very important in complex networks. They 
can not only characterize the detail structure of networks 
but also relate to the structural correlations, motifs, ro- 
bustness and redundancy of pathways, and affect some 
dynamical as well as equilibrium critical phenomena of 
the networks [1,2]. Recently, to avoid the effect of the 
loop structure, some researchers even study the acyclic 
networks [3, 4]. 

For the self-avoiding loop, researchers focus mainly on 
two aspects: the total number of loops and the the dy- 
namic effect of the loop structure. In the former case, 
many counting methods have been proposed [5- 11]. In 
undirected networks, short loops can be exactly counted 
in terms of powers of the adjacency matrix[8]. This 
method can not deal with long loops, because the count- 
ing equation will become very complicated when facing 
long loops. In directed networks, short loops can be es- 
timated by using Nl — jTrA L while long loops can 
be calculated by the entropy[9,10]. Moreover, some re- 
searchers analyze statistics of loops with different length 
L. They use the Monte Carlo sampling method to get 
the frequency and find that the loop number is sharply 
peaked around a characteristic loop length L*. Also, 
they use L* and the relevant index to characterize the 
networks [11]. On the other hand, the dynamic effect of 
the loop structure has been studied frequently. It has 
been pointed out that the loop structure is related to the 
activity in neural networks such as self-sustained activ- 
ities^, 13, 14] and synchronization^, 15, 16]. Specifically, 
the self-sustained activity can not survive without the 
loop structure and the synchronization will be weakened 
when emerging a dominant loop in the network. What's 
more, a scaling behavior of loops is used to explain some 
critical phenomenon in percolation [17] and loop number 
is also used as a ranking method to quantify the role of 
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both nodes and links [18]. 

However, many problems about loops still remain un- 
noticed. In the ref.[10], Bianconi, Gulbahce and Mot- 
ter find that many real networks have fewer loops than 
the counterpart random networks which are a kind of 
random networks with the same number of in- and out- 
links in each node as the real networks. Actually, the 
loops number in different parts of a network varies ac- 
cording to the function of the regions. For example, 
the feed-forward part of the neural networks are sparse 
of loops [3] while other parts in the brain need loops to 
carry out self-sustained oscillation for precessing informa- 
tion^, 14]. For the food web networks, the loops num- 
ber in the metazoan part are relatively small while there 
are many short loops among the microorganisms, called 
microbial loops, for fixed carbon repacking and recovery 
path of ecosystem [19]. Obviously, loops locate unevenly 
in many real networks. Some communities of these net- 
works will be rich in loops while loops will be sparse in 
other parts. This leads us to an interesting question: 
what is the detail organization of loops location like in 
the networks? In this paper, we focus on the unevenness 
of loops location. We first study the distribution of loops 
on single nodes. Then we analyze the rich loop core phe- 
nomenon of uneven loops location by a node removing 
process in some real networks. Finally, we propose an 
index to measure the unevenness. 

Heterogenous distribution of loops on single 
nodes. In the first step, we should study the loops on 
each single node to help us understand how the loops lo- 
cate in the network. For a given network with size N, 
if we want to obtain how many loops passing through a 
specific node, we can simply remove the node from the 
network and count how many loops decreases, the decre- 
ment is the number of loops on this node. As Nl is the 
loops number with the length I of a network, we de- 
note the Nl(i) as the number of loops with length L in 
the network after the node i is removed. So the node 
has Cl{i) = Nl — Nl(i) loops with length L passing 
through. The number of short loops in directed net- 
work can be expressed in terms of powers of the adja- 
cency matrix. In particular, Nl — -^TrA L , provided 
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cause Cl(i) — Nl — iVi(i), we can count the short loops 
on single nodes in any networks. 

Here, we focus on the distribution of Cx(i) in different 
networks. We can compare C£(z) of real networks with 
C£(i) of the counterpart random networks and C£(i) of 
the corresponding ER random networks. The counter- 
part random networks are a kind of uncorrelated random 
networks with the same number of in- and out-links in 
each node as the real networks. The corresponding ER 
random networks is given with the same size and the total 
number of links as the real networks. Of course, C£(i), 
C£(i) and C e L {i) can be obtained by C L {i) = N L -N L (i). 
Actually, the expected value of C£(i) can be gained by 
the formula based on the degree sequence. Motivated 
by the formula of the expected number of loops in the 
uncorrelated random network[10,20], we derive the ex- 
pected number of loops on single nodes in undirected 
and directed random networks. 

For undirected random networks, the expected number 
E(Nl) of short loops with length L is given by [20] 

where k is the degree sequence of the network and < . > 
represents the average value of a sequence. We can obtain 
the expected number of short loops on a specific node as 

E{C L {€)) = ~ A( (B "* < it , -"g" 4fc,) ) L . ( 2 ) 

N N 

where a = Y, kh{kh — 1) and b — ^h- 

h=l h=l 

For directed random networks, the expected number 
E(Nl) of short loops with length L can be obtained 
by[10] 

E{NL)= k { < <i° n i > )L - (3) 

Like the undirected network, we also deduce a formula to 
estimate the expected number of short loops on a specific 
node i. The formula of E(Cl(i)) is 

N N 

where c = k™k^ ut and d = ■ To examine the 

h=l h=l 

validity of our formula, we calculate the exact short loops 
number [8] in directed and undirected random networks 
with prearranged poisson degree sequences and compare 
them with the expected values from our formulas. The 
result shows our formulas can perfectly predict the Cz(i) 
of both undirected and directed random networks with 
given degree sequences. So far, given the degree sequence 
of directed and undirected networks, we can use these two 
formulas to predict the loops number on each node in the 
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FIG. 1: The zipf plots of C r L (i), C£(i), C£(i) and £(C£(i)) 
of four different networks including (a) C.elegans' neural net- 
work, (b)Littlerock food web network, (c)High technology 
company employees' friendship network, (d)Stmark food web 
network. In this Figure, we use L = 5 as example. The links 
for the random networks are averaged by 100 times. 
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FIG. 2: The rank clocks of Cl(i) of C.elegans' neural network, 
with each axis running from rank 1 on the circumference to 
306 at the center. The top rank 1 represents the node with the 
most loops passing through while the bottom rank 306 is the 
node with the fewest loops. Each vertex on the circumference 
stands for the length of the loops. 

uncorrelated random networks. That is to say the dis- 
tribution C£(z) can be represented by the distribution of 
E(C1(i)) which can be simply calculated by our formu- 
las. 

In fig.l, we compare the zipf plots of C£(i), C£(i), 
C£(i) and i?(C£(i)) in four different networks including 
the C.elegans' neural network, the littlerock food web 
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network, high technology company employees' friendship 
network and the stmark food web network. If a node has 
no loop passing through, Cl(i) — 0. In the zipf plot, this 
node will not appear in the log axis. Hence, the shorter 
tail of the line means all the loops are inclined to locate 
in several specific nodes. In addition, the steeper slope of 
the line in zipf plots indicates the distribution of Cl (i) is 
more skewed, which means that nodes are quiet different 
from each other in loops number. If these two features 
are more significant, the distribution of the Cl(%) will be 
more heterogenous. Comparing the real networks to ER 
random networks, we find that the loops are more het- 
erogenous in some real networks such as neural networks 
and some food web networks. Moreover, from the C£(i) 
and C£(i), we can easily find that the degree sequence 
are not sufficient to describe the heterogeneity of loops 
distribution on nodes. In many cases, C£(i) performs a 
more significant heterogeneity than £?(C£(i)) and C£(i). 
However, loops are distributed in some social networks 
almost the same as in the counterpart networks. A typi- 
cal example is given in fig. 1(c). 

In particular, we study Cx(i) with different lengths in 
the C.elegans' neural networks. In fig. 2, we use the rank 
clocks to test whether the rank of loops on each node 
varies dramatically with different loops length £[21]. In 
fig. 2, the rank clocks show perfect pentagon in the cir- 
cumference which means that the top ranks of Ci,(i) do 
not change. Hence, we can clearly see that although for 
different length £, the heterogenous loops distributions 
share the same top rank. 

From the heterogenous distribution of the loops in each 
node, we can easily find that loops locate more uneven 
in some real networks than in counterpart random ones. 
That is to say, some nodes of the real network are rela- 
tively rich in loops while loops are sparse in some other 
nodes. This phenomenon indicates that besides the total 
number, the detail organization of these loops in the real 
networks is quiet different from the counterpart random 
networks. If these nodes with many loops tightly con- 
nect with each other in the same community, of course 
this community will be extremely rich in loops. Com- 
bined with the result in fig. 2, the same top rank in the 
heterogenous loops distribution will enhance the richness 
in loops in the community. In the following section, we 
will discuss the phenomenon of uneven loops location by 
studying the short loops in a specific kind of community 
of networks. 

The rich loop core phenomenon. To investigate 
the detail organization of the self-avoiding loops in a net- 
work, we will study the loops number of different commu- 
nities. In this paper, we consider that a loop belongs to a 
community only if all the nodes of the loops are included 
in this community. For these real networks we are about 
to analyze in this paper, because long loops are strongly 
related to the size of a chosen community, short loops are 
the main elements in the community. Therefore, we only 
consider short loops in this paper. It has been investi- 
gated that many real networks have fewer loops than the 




FIG. 3: The change of the loops number as the nodes remov- 
ing from three kinds of networks. In this figure, we choose 
L = 5 as a example. In fact, because the community ob- 
tained in this process is the selected community, the loops 
number in these communities is Ml(ti). Furthermore, the 
size of selected community n = N — t in each step t. The 
links for the random networks are averaged by 50 times. 
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FIG. 4: The loop density ratio J/z,(n) of three network in each 
step of node removing process. We choose L = 5 as a exam- 
ple. It is clear that loops locate more unevenly in C.elegans' 
neural network than the counterpart random networks. Fur- 
thermore, the size of selected community n — N — t in each 
step t. The links for the random networks are averaged by 50 
times. 

randomized counterparts in both short loops and long 
loops such as C.elegans' neural network, Food Web net- 
works, Power-grid networks etc[10]. Here, we take some 
of these real networks to make a further study of the 
loops in their communities. 

First, we introduce a node removing process. In each 
step t of the node removing process, we remove the node 
with the smallest Cx(z) hi the network. It means that the 
node with fewer loops passing through will be removed 
first. Specifically, the Cx(j) should be updated in each 
step. After N steps, we can remove all the nodes from 
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FIG. 5: Typical results of comparing selected communities' 
loops of many networks to that of the counterpart random 
models. The results are averaged by 50 times, (a) C.elegans' 
neural network, (b) C.elegans' metabolic network, (c) and (e) 
are food web networks in two different places, (d) Company 
employees' friendship network which is obtained by question- 
naire, (f) directed ER random networks with 100 nodes and 
800 links. 



the network. If we want to obtain a community with 
size n, the number of the removed nodes should be t = 
N — n. Because the node about to be removed in every 
step has the fewest loops among the remaining nodes, 
the community obtained by this process has the richest 
short loops among all the communities with the same 
size. The community obtained by the node removing 
process is called the selected community in this paper, 
the loops number in this community is denoted as Ml in) . 
Obviously, M L {N) = N L and M L (0) = 0. If a network is 
greatly uneven in the loop location, Ml(ti) will decrease 
slowly in the beginning while decline dramatically in the 
end through out the node removing process. 

We analyze the C.elegans' neural network by the node 
removing process as an example. In fig. 3, we compare 
the C.elegans' neural network, the counterpart random 
networks and the corresponding ER random networks. 
It can be seen that although the counterpart random 
networks have more short loops than C.elegans' neural 
network, it can not compete with the C.elegans' neu- 
ral network in some specific selected communities. For 
instance, the selected community with 50 nodes in the 
C.elegans' neural network has more short loops than that 
in the counterpart random networks as in fig. 3. 

Furthermore, we define the loop density as the propor- 
tion of the loops number and the community size. So we 
can easily get the loop density of the original network 
Nl/N. To compare the loop density of the selected com- 
munities and the original network, we define the loop 



density ratio rji J (n) in each selected community during 
the node removing process as 



M L {n)N 



(5) 



The loop density ratio r\L (n) is the proportion of the 
selected community's loop density and the original net- 
work's loop density. If the loops locate unevenly in a 
network, T]L[n) will become larger than 1 during the 
node removing process. It indicates some communities 
in the network have bigger loop density than the original 
network. Obviously, in the C.elegans' neural network, 
some of the selected communities enjoy larger loop den- 
sity than the original network, so rii,(n) is larger than 
1 as shown in Fig. 4. Although the counterpart random 
networks can also have a r)L(n) larger than 1, its value 
is always lower than that of C.elegans' neural network. 
This feature indicates that loops locate more unevenly in 
C.elegans' neural network than in the counterpart ran- 
dom networks, which means that lots of loops are limited 
in some specific nodes of the C.elegans' neural network. 

So the unevenness may result in a rich loop core phe- 
nomenon, which means some selected communities with 
far higher loop density compared with the original net- 
works. Although many real networks have fewer short 
loops compared to counterpart random models, the rich 
loop phenomenon will make some communities in real 
networks more loopy than the corresponding communi- 
ties in the counterpart random networks. It can be de- 
tected by the ratio as: 



Pl{u) 



M real 



(n) 



M r L and 
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where n is the size of the selected community. For ex- 
ample, the index will turn from pl(it-) < 1 to Ph{n) > 1 
as nodes removed in the C.elegans' neural network as 
show in fig. 3. Also, we also investigate many other real 
networks, some typical results are shown in fig. 4. 

As in fig. 4 (a), (c) and (e), although these real net- 
works, such as neural networks and some food web net- 
works, have fewer total loops than the counterpart ran- 
dom networks, some selected communities of them are 
more loopy than that of counterpart random ones. More- 
over, we find that short loops with different lengths L per- 
form the same trend in the rich loop core phenomenon. 
If the loops with specific length L locate unevenly in 
the networks, loops with other lengths locate uneven as 
well. However, not all the networks have this kind of phe- 
nomenon as shown in Fig. 4(b). Some real networks with 
far fewer loops than the counterpart random networks 
have pL(n) < 1 for all n. For example, the C.elegans' 
metabolic network belongs to this category. In addition, 
we find that some social networks have more short loops 
than the counterpart random ones as in fig. 4(d). These 
networks are not discussed in ref[10]. This category in- 
cludes the prisoners' friendship network, high technology 
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employees' friendship network, the family visit network, 
the flying-team partner choosing network, the dining ta- 
ble partner choosing network and so on [22]. If these net- 
works enjoy more uneven loops location, their selected 
communities can only be more and more loopy than that 
of counterpart random networks as the community size 
n varies, see fig. 4(d). Finally, we use the index Ph(n) 
to detect a ER random network with 100 nodes and 800 
links. The result in fig. 4(f) shows that Pl(ti) = 1 ap- 
proximately. 

Besides the rich loop core phenomenon, how to find the 
rich loop core in a network is an interesting question. In 
this paper, we simply consider the rich loop core appears 
at the maximum t]l (ni) j which means the rich loop core 
will have highest loop density than any other community 
in the networks. Some typical results are shown in fig. 6. 
The littlerock food web and the C.elegans' metabolic net- 
work have significant rich loop cores which indicates they 
have much higher loop density community compared with 
the original networks. The C.elegans' neural networks 
and the StMark food web have such rich loop cores as 
well. On the contrary, this phenomenon is not so obvi- 
ous in company employees networks and large degree ER 
random networks. For the C.elegans' metabolic network, 
the high loop density ratio is due to the small number 
of total loops. These loops unavoidably locate in sev- 
eral specific nodes, so the loops density in the selected 
community will be very large compared with the original 
network. Furthermore, we analyze the rich loop cores of 
the food web and neural networks. It is interesting that 
most of the nodes in the cores are from interneurons in 
C.elegans' neural network and from microorganisms in 
food web networks. Specifically, the rich loop cores of 
C.elegans' neural networks have different sizes n for dif- 
ferent loop lengths L. For example, n = 31 when L = 3, 
n = 38 when L = 4, n — 50 when L = 5 and n = 51 when 
L = 6. These rich loop cores share 30 nodes and 24 of 
them are interneurons. The C.elegans' neural networks 
are composed by sensory neurons, interneurons and mo- 
tor neurons. Most interneurons are in the nerve ring gan- 
glia. Their main function is to process signals [23, 24]. So 
the loops number in these interneurons is relatively larger 
than the others. Likewise, the rich loop cores in stmark 
food web networks are n — 20 when L = 3, n = 22 
when L = 4, n = 24 when L = 5 and n — 24 when 
L = 6. These rich loop cores share 14 nodes and 10 of 
them are microorganisms. In the food web networks, the 
microorganisms are the main element in microbial loops 
which are strongly related to fixed carbon repacking and 
recovery path of ecosystem[19]. 

Measurement for unevenness of loop location. 

In order to quantify how unevenly loops locate in the 
networks, we present an index which bases on the node 
removing process. In order to Simplify the computing 
complexity, we do not update Cl(i) in each step dur- 
ing the node removing process in this section. We test 
and find the result obtained in way is sufficient to rep- 
resent that by updating Cl(i) in each step statistically. 
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FIG. 6: The rich loop core phenomenon of some typical real 
networks. The links for the random networks are averaged 
by 50 times, (a) C.elegans' neural network, (b) C.elegans' 
metabolic network, (c) and (e) are food web networks in two 
different places, (d) Company employees' friendship network 
which is obtained by questionnaire, (f) directed ER random 
networks with 100 nodes and 800 links. 

Moreover, we remove the nodes based on the attacking 
rate p. For example, if p = 0, no node is removed and 
the network is all the same with the original network 
and the loop number is Nl. Hp — 0.1, we just remove 
[pN] ce n nodes from the network and the loop number is 
Ml(N —\pN] ce n). Here, [] ce u represents the operation of 
rounding upward. Then, we use Al(jp) — MldJlzMlksill 
to normalize the loop number of each community so that 
Al(j?) which is corresponding to the loop number declines 
from 1 to through out the node removing process as 
shown in fig. 7. Again, we use L = 5 as an example in 
ng.7. 

If a network is significantly uneven in the loop location, 
Al(p) will decline slowly in the beginning while dramat- 
ically in the end during the node removing process. On 
the contrary, if loops locate evenly in the network, the 
Al(p) will decline almost the same as the correspond- 
ing ER random networks. Therefore, the unevenness of 
loop location can be measured by the difference between 
the real network and the corresponding ER random net- 
work. Here, we use /£(p) = A r L (p) — A e L (p) to estimate 
the difference. So the unevenness of loops location can 
be represented by the index as 

Rl= f I L (p)d P , (7) 
Jo 

where — 1 < Rl < 1. The severer the unevenness is, the 
larger the index Rl is, which means the loop location 
departs more largely from the corresponding ER random 
network. Of course, the index Rl can also be used in 
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P P 



FIG. 7: the A r L (p), A c L (p) and A L (p) change as the p in differ- 
ent networks including C.elegans' neural network, high tech- 
nology company employees' friendship network, littlerock FW 
network, and ER random network with small degree. L is 
chosen as 5 and the result is averaged by 50 times. 



analyzing the counterpart random network. 

Actually, the index Rl is the area between the lines of 
the real network and the ER random network in hg.7. It 
can be seen that the real networks and the counterpart 
random networks can be different in the unevenness, as 
the index R L ^ R C L . Typically, if the real network is 
very sparse in the links, the network has only small num- 
ber of loops. This will lead to a phenomenon that some 
part of the network has some loops, while other part has 
no loop at all. However, the corresponding ER random 
network has the same condition too. For the Rl is the 
area between the lines of the real network and the ER 
random network, the R^ will be a small value under this 
circumstance. Whether this uneven location of loops re- 
sults from the specific structure of the real network or 
from the small degree as the ER random network can be 
estimated from the index R^. As the figure(d) in the 
fig. 7, we use the index Rl to investigate a ER random 
network with small degree as K = 2. The result shows 
that the three lines are almost the same meaning that 
R L « 0. 

Additionally, we consider several more directed real 
networks[22]. As mentioned above, although each net- 
work has different kinds of short loops based on length 
L, these loops perform almost the same. We use aver- 
age Rl to represent the unevenness in loops location. It 
can be gained by R =< Rl > and the L = 3, 4, L max 
where L max — 8 according to Ref[9]. The index R for 
these real networks are given in table 1. From the R r , 
how significant the unevenness in loops location is can 
be known. By comparing the R r and the R c , we can 
distinguish whether this unevenness results from the de- 
gree sequence. In table 1, it can be seen from the R r 
that the C.elegans neural network, C.elegans metabolic 



TABLE I: Results of the analysis of networks based on index 
R 



network 


size 


links 




Rc 


Rr — Rc 


C.elegans' neural 


306 


2359 


0.436 


n 


31 / 


0.119 


C.elegans' metabolic 


453 


2040 


0.468 


n 
U 


OO ( 


0.131 


E.coli's metabolic 


896 


958 


-0.021 




U.UUZ 


-0.019 


Mondego FW 


46 


400 


0.231 





257 


-0.026 


Michigan FW 


39 


221 


0.194 





207 


-0.013 


Littlerock FW 


183 


2494 


0.672 





439 


0.233 


StMarks FW 


54 


356 


0.378 





256 


0.122 


Prisoners 


67 


182 


0.169 





004 


0.165 


Flying-teamers 


48 


351 


0.108 





060 


0.048 


Company employees 


36 


147 


0.219 





163 


0.056 


ER random 


















network and some food web networks are really uneven 
in loop location. The social networks, Escherichia Coli's 
metabolic network and some other food web networks do 
not have such significant unevenness. In fact, the number 
of different species will affect the loop location in food 
web networks. For example, too many microorganisms 
will make the loops more even and too many metazoans 
will reduce the total loops number, at this time the rich 
loop core will be more obvious. Additionally, both the 
C.elegans' metabolic network and the Escherichia Coli's 
metabolic network have very few loops, but the degree 
sequence and the total links of the former one allows 
the counterpart random and the ER random networks 
to have much more loops while the latter one does not. 
Hence, the R r of these two metabolic networks are dif- 
ferent. Moreover, comparing the R r and the R c in table 
1, it can be found that the degree sequence is not suffi- 
cient to describe the unevenness in loops location. It is 
clear that the C.elegans' neural network, littlerock and 
stmark food web networks are more uneven in loops lo- 
cation than the counterpart networks. That is why they 
can have some communities more loopy than the corre- 
sponding communities in the counterpart random models 
despite the fewer total loops number. 

Conclusion. The previous works on the loops mainly 
focus on the total number of loops and the dynamic effect 
of the loop structure. However, the loop location is also 
very important in networks. Generally, loops tend to 
locate in some specific nodes in some real networks, which 
means some communities of the network are extremely 
rich in loops while the loops are relatively sparse in other 
parts. If this uneven location is significant enough, the 
rich loop core phenomenon can be formed in some real 
networks. 

The rich loop core phenomenon is meaningful for the 
typical function of real networks. For instance, the loop 
structure is strongly related to the self-sustained activ- 
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ities in neural networks, so the rich loop core may help 
to understand the functional regions in the neural net- 
works. For the food web networks, almost all the loops in 
rich loop cores of food web networks are microbial loops 
and plays an important part in fixed carbon repacking 
and recovery path of ecosystem. In addition, this uneven 
location of loops may provide a new way to study the 



community detection in directed networks, which asks 
for further research. 
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